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Abstract:-In the previous paper, the use of matrix iteration to determine the natural frequencies of vibration of 
continuous beam system using the concept of wave propagation in a prismatic bar is reported. The natural 
frequencies of vibration obtained from the formulated model were compared with the values obtained from 
literature and those of the exact solution. An error of 28% and 30% was observed when the fundamental frequency 
of vibration obtained from the previous model was compared with that obtained from literature and the exact 
solution. In the present study, an improved mathematical model is formulated to determine the fundamental 
frequency of vibration of a continuous beam as a structural system with distributed mass using mode shape 
concept. In the course of the system's oscillation, the displacements produced by the force of inertia is assumed to 
have the shape of a particular vibration mode and in harmonic with the particular modal frequency. A numerical 
example was given to demonstrate the applicability of the present model. The fundamental frequency value 
obtained from the present model was found to improve by 27.4% and almost identical to the exact fundamental 
frequency value and that obtained literature. 
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I. INTRODUCTION 

Resonance is a very dangerous phenomenon that occurs if the fundamental value of the natural frequencies of 
vibration is exceeded by the frequency of excitation. Resonance results in large amplitude displacement of structures leading 
to development of large stresses and strains in the affected structure that may eventually result to structural failure thereby 
affecting the performance of the structure or structures in service. The dynamic analysis of a continuous beam as a system 
with distributed mass to get the fundamental vibration frequency is always cumbersome because of the difficulty involved in 
the mathematical manipulations. The difficulty that is encountered in the mathematical process is due to infinite number of 
degrees of freedom [1-9]. The dynamic analysis is made simpler using lumped mass concept. The original beam with 
distributed mass is now converted to a weightless system with masses lumped at chosen points called the nodal points. The 
weightless beam system now has a finite number of degrees of freedom [10]. The degree of freedom is numerically equal to 
the number of independent geometric parameters that describes the positions of all masses for all possible displacements of 
the structural system at any point in time. The present model is said to be defined if the nodal lumped masses and their 
coordinates are known [11]. 

In this paper, matrix iteration is employed to determine the fundamental frequency of vibration of a continuous 
beam system undergoing self excited vibration. The algorithm involved is simple and can be achieved manually most 
especially when finite number of degrees of freedom is involved. 

II. 2. Formulation of Mathematical Model 

For an undamped MDOF beam system (Figure 1) with n degrees of freedom, displacements are assumed to be 
linear functions of the forces of inertia, the equations of motion of the lumped masses at chosen nodal points are given by: 

- X x =m l X 1 S n +m 2 X 2 S u + m n X n 8 u 



~ X 2 = m 2 X 2^21 +^2^2^22 + m „ X n S 2n 



-X„ =m„XS. +m X S„ + m„XS„„ 

n n n nY I I nl n n nn 



(i) 



where: 

mX . = Inertia force generated by an ith particular oscillating mass. 

X t = displacement produced at ith mode by inertia force generated by an ith oscillating mass. 

Of, , O u = direct and indirect flexibility coefficient respectively. 
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The flexibility coefficients resulting from the forces of inertia at the individual nodal points are given by: 



n — ir i u 

o mi m 2 I, mi ^ m N „ 

Figure 1: Lumped masses at beam nodes. 

^ ;; ='V — [ m.m.dx (2) 

The equation for undamped natural vibration frequency is given by : 

K = co 2 m (3) 

For an ith vibration mode equation (3) transforms to: 

K = cofm (4) 

Multiplying both sides of equation (4) by an ith displacement X; caused by an ith force of inertia X . transforms equation 

(4) to: 

KX^cofX.m (5) 

The displacement that is produced by an ith oscillating mass has the configuration of an ith vibration mode and also in 
harmonic with an ith modal frequency. Therefore, 

KX t =F t (6) 

Equation (5) now becomes: 

F^afX.m (7) 

From equation (6), 

K 

where: 

F, = Force of inertia at ith mode 

K Inverse of stiffness matrix 



X^^FK- 1 (8) 



Substituting for F t in equation (8) using equation (7) gives: 

X. =a)?mK- l X i (9) 



Let the row vector of n dimensions 

x t =[x lt x 2t ... t xj ( 10 > 

represent the vector of displacement in the ith vibration mode. 
From structural mechanics, 

K l =8 (ii) 

where: 

O = flexibility coefficient matrix. 

Without loss of generality, 

K? = 8 tt (12) 

and 

K g l =S, (13) 

Equation (7) now transforms to: 

X t -cofmSX i (14) 

Let 

mS =h (15) 

Equation (13) represents the dynamic charcateristics of the weightless beam subjected to oscillating masses. 
Therefore, from equation (12), 

X. =a>fhX i (16) 

For the first vibration mode i=l, and equation (5) becomes: 
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X x =a\hX° (17) 

where: 

x i ; =a 1 if (is) 

represents an arbitrarily chosen initial displacement vector. 

Let y represent non-zero nth order displacement vector and let max y max represent its largest displacement element. Let 

X be the vector obtained by when the entries of y are scaled by y max . 
Using equation (16), 

y„ 



y 
The generalized sequence of improved displacement vectors are given by equations (18) and (19). 

y k 



x 



•^^ma; 



-; A; = 0,1, ...,n 



hX k =h T X k =[h] 



A\ 



x„ 



: y k , 



Let i^t/i t> e a se qu ence °f approximations to y n with 

lim (y k )=K 

k=> oo 

where: 



\ = co 2 n (n = 1,2,-) 

From equations 
®„=fat + l)°-' 



From equations (19) and (21), the natural frequency corresponding to a given vibration mode is given by: 
).5 



and the fundamental value (n=l) is given by 

»l=U + l) ' 5 

From statical consideration, the ith modal mass at ith nodal point over an ith weightless beam segment is given by: 



m, 



where: 



P 



fe+0 



(19) 



(20) 



(21) 



(22) 

(23) 
(24) 
(25) 

(26) 



p = distributed mass intensity of the beam in Kg/m. 

III. CONCLUSION AND RESULTS 

An Example for Numerical Study 

A numerical example is used to demonstrate the applicability of the present formulation. A simple supported 
uniform beam having a distributed mass intensity of 4.75Kg/m as shown in Figure 2 is used for this numerical study[l]. 



mi 
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Figure 2: A 3 degrees of freedom beam system for numerical study. 
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P = \ 

m, rn 2 m 3 





Figure 2: Derivation of flexibility factors 



The flexibility matrix is symmetric. Therefore, 

8.. =S.. 

U J' 

Using equation (2), the flexibility factors are obtained as follows: 
rl 2 

EIS n =1 m l dx 
_ 1.070 



EI 



ES 22 =\ m\dx 
2.604 



S. 



22 



EI 



1 rl 2 

S 3J = — m 3 dx 



Ei- 




4.0m 



-y 



0.8 + 0.8 



X- 



2.5m 



-y 




1.25+1.25 



0.8 + 0.8 
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EL 



5 xi = 5 ii = ir T \ m \ m 2 dx '- 

' ' Jo 




1.5m . 



0.8 0.8 



0.5 0.2 




1.480 
EI 




1.5m 

I 1 



0.5 0.5 



1.25 O.f 




s 13 = s 3l 



J_f' 

EI 



\ m l m J - 

Jo 

1.0m 




0.8 0.8 0-5 0.5 



1.0m 




1.5m 

I 1 



1.0m 

I 1 



0.2 0.2 




1.5m 1.5m 

I 1 I 1 



1.0m 

I 1 



0.707 
EI 




0.2 o.: 



0.5 0.5 



0.8 0.S 




1 r' 



^23 = ^32 = — ~\ m 2 m 3 dx= 

HI Jo 



EI- 



2.5m 




h 



1.5m 1.0m 

1 I 1 



1.25 1.25 



0.5 0.5 




1.5m 



1.0m 
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From equation (26), 



34 



Application of Matrix Iteration for Determining the Fundamental. 



m, 



m, 



m. 



p + l.5p 



-1.25 p 



1.5/7 + 1.5/7 



= 1.5/7 



The flexibility coefficients at the three nodal points are now arranged in matrix form as follows: 

1.070 1.480 0.770" 

1.480 2.604 1.480 

0.770 1.480 1.070 

From equation (15), the dynamic matrix is given by: 



J_ 

EI 



1.070 


1.480 


0.770 


1.25/7 





1.480 


2.604 


1.480 





15/7 


0.770 


1.480 


1.070 











1.25/7 



1 
EI 



'1.34/? 2.22/7 0.96/7 
1.85/7 3.91/7 1.85/7 
0.96/7 2.22/7 1.34/7 



Using equation (18), multiplication of dynamic matrix with assumed initial unit displacement vector gives: 



X, 



1 

EI 



1.34/7 2.22/7 0.96/7 
1.85/7 3.91/7 1.85/7 
0.96/7 2.22/7 1.34/7 



1 
£7 



4.52/7 
7.61/7 
4.52/7 



Using 7.61/7 as the maximum displacement in the non zero displacement vector, the improved displacement vector is 
given by: 

X 2 =— (3.59/7,6.10/7,3.59/7) r 
EI 

Again, the maximum displacement is 6. 10/7 . The new improved displacement is: 

X 3 =— (3.575/?, 6.09/7, 3.575/?) r 
EI 

Using equation (23), 

[ei 



CO, 



\-*L- =0.405 
6.09/7 \ p 



Fundamental 
frequency 



Table 1: Comparison of results 



Present model 



0.4050 



Sule 2009 



0.515 



Osadebe 1999 



0.4039 



Exact solution 



0.3948 



IV. DISCUSSION OF RESULTS 

Table 1 shows the comparison of result obtained from the present model with Sule [2], Osadebe [1] and the exact 
solution [3]. The percentage errors of 27.5% and 30.45% in the previous model[2] compared with the control points [1] and 
exact solution [4] have reduced to 0.272% and 0.258% in the present formulation showing the effectiveness of the present 
model in the determination of the fundamental frequency of vibration of a continuous beam. The disparity between the result 
of the present model and the previous model[l], Osadebe [2] and the exact solution[3] may be due to difference in the 
assumptions used in the model formulation. 

V. CONCLUSION 

In conclusion, the present model produces a result that is almost identical with those of the control points [1] and 
[4] and improves on the previous result of fundamental frequency [2] by 27.4%, showing higher predictive ability of the 
present model. The present model can be used to predict the fundamental frequency of vibration of a multi-storey building. 
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